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SECTION  I 


INTRODUCTION 


1.1  BACKGROUND. 

The  use  of  impedance  testing  and  its  applica¬ 
tion  to  buried  arch  response  is  discussed  in  this 
report.  A  goal  of  impedance  testing  is  to 
produce  functional  relations  which  character¬ 
ize  the  shock  input  levels  at  points  of  equipment 
attachment  to  the  structure.  The  experimentally 
determined  impedances  of  the  structure  are 
used  to  develop  a  transfer  function  between  the 
generalized  input  loads  upon  the  structure  and 
the  dynamic  response  velocities  at  the  attach¬ 
ment  points  of  equipment.  The  restriction  of 
this  methodology  is  that  linear  behavior  of  both 
the  structure  and  the  soil  must  exist.  The  reality 
is  that  some  degree  of  nonlinearity  occurs  at 
any  significant  pressure  level  with  the  onset  of 
soil  nonlinearity  preceding  the  onset  of  struc¬ 
tural  nonlinearity. 

The  particular  impedance  tests  studied  in  this 
report  are  the  forced  vibration  tests  performed 
on  uncovered  and  shallow  buried  arch  and 
rectangular  structures  in  the  ESSEX  V  Phase  3 
test  series.  (Ref.  1) 

For  these  tests,  the  original  backfill  was  exca¬ 
vated  and  the  uncovered  structure  was  sub¬ 
jected  to  forced  vibration  tests.  The  sand 
backfill  was  then  replaced  by  raining  it  from  a 
height  of  five  feet  and  the  forced  vibration  tests 
were  repeated.  Figure  1  shows  the  geometry  of 
the  arch  structure  and  the  backfill  and  instru¬ 
mentation  layout. 

Test  results  include  impedance  (Fmax/Vmax) 
versus  frequency  relationships.  Resonance  is 
defined  to  occur  when  the  driving  frequency 
approaches  a  particular  natural  frequency  of  the 
structure.  This  is  exhibited  by  a  valley  in  the 
impedance  versus  frequency  curve. 


The  test  results  showed  that  for  the  rectangular 
structures,  resonance  was  observed  in  the 
shallow  buried  structures  at  the  same  frequen¬ 
cies  as  in  the  uncovered  box  structure.  For  the 
arch  structure,  resonance  was  not  experimen¬ 
tally  observed  in  the  shallow  buried  arches. 

In  an  attempt  to  explain  the  differences  in  the 
experimental  results  of  the  boxes  and  the 
arches,  linear  time  domain  analyses  were 
previously  performed  to  simulate  the  behavior 
of  these  structures  (Ref.  2).  These  analyses 
examined  what  levels  of  radiation  damping 
would  be  needed  to  supress  resonance.  It  was 
determined  that  for  the  box  structure,  radiation 
damping  of  the  order  of  10%  to  20%  was 
needed  to  supress  resonance.  For  the  arch 
structure,  it  was  determined  that  a  radiation 
damping  level  of  50%  would  be  needed  to 
supress  damping.  Since  measured  values  of 
radiation  damping  are  in  the  range  of  10  to  1 5%. 
it  was  concluded  that  radiation  damping  alone 
could  not  account  for  the  supression  of  damp¬ 
ing  for  the  covered  arches. 

Identified  in  Ref.  2  as  an  additional  damping 
source  was  hysteretic  behavior  in  the  soil 
adjacent  to  the  structure.  Hysteresis  in  soil 
stress-strain  relations  describes  the  energy 
losses  produced  by  the  load/unload  and  reload 
cycles.  These  cycles  are  produced  in  the  soil  by 
the  driving  oscillations  of  the  buried  structure. 
The  contributing  effect  of  the  soil  hysteresis 
upon  the  experimentally  observed  supression 
of  arch  resonance  will  be  further  discussed  in 
this  report. 

1.2  SCOPE  AND  OBJECTIVES. 

This  report  discusses  the  development  and 
validation  of  analytical  techniques  for  hardness 
assessment  of  shallow  buried  arches.  It  repre¬ 
sents  an  extension  of  previous  studies  done  by 


1 


Weicllinger  Associates  under  DNA  Sponsor¬ 
ship  to  simulate  the  forced  vibration  testing  of 
uncovered  and  covered  arches  and  boxes  in 
Project  ESSEX.  Those  studies  indicated  that 
material  and  geometrical  non-linearities  need¬ 
ed  to  be  considered  when  modeling  shallow 
buried  arches. 

The  role  of  the  various  component  soil  material 
non-linearities  and  their  interplay  with  the  arch 
structure  is  examined  in  this  report.  Particular 
attention  is  drawn  to  the  responsibility  of  soil 
non-linearities  in  helping  to  describe  the  appre¬ 
ciable  difference  between  the  shallow  buried 
arch  and  box  behavior. 

Time  marching  simulation  using  non-linear 
finite  element  analysis  is  performed  to  over¬ 
come  the  theoretical  obstacle  presented  by  the 
non-linear  soil-structure  system.  Various  ana¬ 
lytical  models  of  the  soil  media  are  used  to 
isolate  the  influences  of  soil  behavior  upon  the 
response  of  buried  arches  to  impedance  testing. 

The  examination  of  soil  response  is  within  the 
scope  of  the  study  of  impedance  testing  as 
applied  to  buried  arches.  Determining  the 
extent  of  non-linear  soil  response  will  greatly 
aid  the  investigation  to  determine  the  extent 
non-linearities  render  the  impedance  technique 
impractical  for  shallow  buried  arches.  This 
study  will  also  assist  in  the  development  of 
analytical  techniques  which  utilize  the  results 
of  impedance  testing  to  develop  shock  input 
levels  at  equipment  attachment  points. 

1.3  IMPEDANCE  TESTING:  METHOD¬ 
OLOGY  AND  RESULTS. 

The  use  of  impedance  testing  and  its  applica¬ 
tion  to  buried  arch  response  is  discussed  in  this 
report.  A  goal  of  impedance  testing  is  to 
produce  functional  relations  which  character¬ 
ize  the  shock  input  levels  at  points  of  equipment 
attachment  to  the  structure.  The  experimentally 
determined  impedances  of  the  structure  are 
used  to  develop  a  transfer  function  between  the 


generalized  input  loads  upon  the  structure  and 
the  dynamic  response  velocities  at  the  attach¬ 
ment  points  of  equipment.  The  restriction  on 
this  methodology  is  that  linear  behavior  of  both 
the  structure  and  the  soil  must  exist.  The  reality 
is  that  some  degree  of  nonlinearity  occurs  at 
any  significant  pressure  level  with  the  onset  of 
soil  nonlinearity  preceding  the  onset  of  struc¬ 
tural  nonlinearity. 

Impedances  are  measured  by  applying  an  input 
forcing  function  to  the  structure  using  either 
electromagnetic  (low  force  and  high  frequen¬ 
cy)  or  electro  hydraulic  (high  force,  low 
frequency)  vibrators.  A  variable  sine  sweep 
rate  is  used  when  driving  the  structure.  The 
impedance  is  determined  using  a  measurement 
point  and  a  series  of  drive  points  on  the 
structure.  The  response  of  the  measurement 
point  to  different  driving  frequencies  is  found 
by  summing  the  products  of  experimentally 
determined  impedance  functions  and  drive 
functions  in  the  frequency  domain. 

The  particular  impedance  tests  studied  in  this 
report  are  the  forced  vibration  tests  performed 
on  uncovered  and  shallow  buried  arch  and 
rectangular  structures  in  the  ESSEX  V  Phase  3 
test  series.  Constant  force  sinusoidal  frequency 
sweeps  were  ran  at  three  different  force  levels. 
Mechanical  impedance  (Fmax/Vmax)  versus 
frequency  plots  were  made.  Resonant  frequen¬ 
cies  were  determined  from  such  system  curves. 
An  outline  of  the  test  plan  for  ARCH  3A 
(Covered  and  Uncovered)  is  below. 

The  test  results  obtained  showed  that  for  the 
rectangular  structures,  resonance  was  observed 
in  the  shallow  buried  structures  at  the  same 
frequencies  as  in  the  uncovered  box  structures. 
For  the  arch  structures,  resonance  was  not 
experimentally  observed  in  the  shallow  buried 
arch  structures.  In  an  attempt  to  explain  the 
differences  in  the  experimental  results  of  the 
boxes  and  the  arches,  analyses  were  performed 
to  simulate  the  behavior  of  certain  types  of 
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shallow  buried  structures  subjected  to  forced 
vibration  testing.  (Refs.  2,3  and  4). 


TEST  DESCRIPTION 

DRIVE  POINT 

FORCING  FUNCTION 

SINGLE  VIBRATOR 

FI 

FREQUENCY  SWEEP 

DUAL  VIBRATORS 

IN  PHASE 

F2  and  F3 

FREQUENCY  SWEEP 
IN  PHASE 

DUAL  VIBRATORS 

OUT  OF  PHASE 

F2  and  F3 

FREQUENCY  SWEEP 
OUT  OF  PHASE 
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SECTION  2 


SIMULATION  OF  BURIED  ARCH  STRUCTURES 


2.1  FINITE  ELEMENT 
METHODOLOGY. 

The  finite  element  method  together  with  an 
explicit  central  difference  integration  scheme 
is  utilized  in  the  analysis  of  shallow  buried 
arches  subjected  to  dynamic  loading.  Both 
material  and  geometrical  nonlinearities  are 
considered  in  the  analysis. 

The  plane  strain  finite  element  program  SAD¬ 
NESS  (SEISMIC  ANALYSIS  for  DYNAMIC 
NONLINEAR  SLOPE  STABILITY)  (Ref.  5) 
is  used.  The  SADNESS  program  was  devel¬ 
oped  to  examine  the  earthquake  induced  failure 
of  earth  slopes.  Its  range  of  applicability 
includes  the  non-linear  response  of  buried 
structures  subjected  to  vibrational  loadings. 
Both  geometrical  and  material  nonlinearities 
are  accounted  for.  Geometrical  nonlinearities 
are  incorporated  through  the  use  of  an  incre¬ 
mentally  updated  Lagrangian  approach  in  con¬ 
junction  with  a  constitutive  equation  relating 
the  co-rotational  stress  rate  to  the  rate  of 
deformation  tensor.  The  theory  is  valid  for  large 
deformations  (strains)  as  well  as  for  large 
deflections  and  rotations.  Material  nonlineari¬ 
ties  are  incorporated  through  the  use  of  a  cap 
model  constitutive  model  extended  to  improve 
its  capability  to  represent  cyclic  behavior 
characteristic  of  sinusoidal  ground  motions. 

The  SADNESS  model  of  the  buried  arch 
(Figure  2)  consisted  of  455  elements  and  267 
nodal  points.  Constant  strain  triangles  are  used 
throughout.  Five  (5)  different  materials  were 
modelled  including  the  concrete  arch,  the 
backfill  and  three  subsequent  soil  layers.  The 
finite  element  mesh  of  the  concrete  arch  is 
composed  of  64  elements.  This  discretization  is 


adequate  to  represent  the  arches  compressive 
circumferential  behavior. 

The  finite  element  mesh  is  extended  in  space  so 
as  to  maximize  the  simulation  time  before  any 
reflected  waves  from  the  finite  element  bound¬ 
aries  reach  the  shallow  buried  arch.  Dashpots 
are  placed  at  the  boundaries  to  absorb  any 
waves  which  are  generated  there. 

For  each  test  of  the  series,  forced  loading  of  the 
form  P  =  Pmax  sin(cot)  is  applied  at  the  drive 
point. 

Direct  integration  in  time  of  the  equations  of 
motion  was  performed  using  a  time  step 
interval  of  60  microseconds.  A  total  of  0.24 
seconds  of  real  time  response  was  simulated. 

The  use  of  the  advanced  soil  plasticity  CAP 
model,  which  satisfies  the  requirements  of 
continuity  and  uniqueness  of  solution  and  has 
the  capability  of  simulating  a  wide  range  of  real 
effects  in  soils  was  deemed  essential  for  these 
analyses  (Refs.  6,  7  and  8). 

The  CAP  model  is  based  on  the  classical 
incremental  plasticity  theory  and  is  capable  of 
representing  the  mechanical  behavior  of  soils 
and  rocks  while  satisfying  all  theoretical  re¬ 
quirements  for  properly  posed  initial/boundary 
value  dynamic  problems.  The  proper  use  of  the 
CAP  model  assures  existence  and  uniqueness 
of  solution  as  well  as  continuous  dependence 
on  the  initial  and  boundary  data.  The  CAP 
model  has  a  non-softening  convex  yield  sur¬ 
face.  For  stress  points  on  the  yield  surface,  the 
plastic  strain  rate  vector  is  outwardly  directed 
and  normal  to  the  yield  surface  in  stress  space. 
The  CAP  model  permits  inelastic  hardening  in 
hydrostatic  loading,  limits  the  amount  of  dilat- 
ancy  during  shear  failure  and  provides  a  good 
fit  to  material  property  data. 
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To  examine  the  interplay  of  the  various  compo¬ 
nents  of  the  soil  material  nonlinearities,  two 
different  soil  CAP  models  are  used.  They  are; 

1.  Elasto-Plastic  soil  behavior  (CAP  Mod¬ 
el) 

2.  Viscoelasto-Plastic  soil  behavior  (Vis¬ 
coelastic  CAP) 

These  models  will  serve  to  isolate  the  role  of 
soil  non-linearities  in  determining  the  response 
of  the  arch  structure  in  forced  vibration  testing. 

2.2  VISCOELASTIC  CAP  MODEL. 

Certain  amendments  to  the  CAP  model  were 
made  to  improve  its  capability  to  represent 
cyclic  behavior  characteristic  of  sinusoidal 
ground  motions.  Much  of  the  previous  use  of 
the  CAP  model  has  been  for  ground  shock 
calculations  involving  much  higher  stress  lev¬ 
els  than  are  applicable  to  low  amplitude 
sinusoidally  induced  motions.  Additionally, 
blast-induced  loadings  generally  include  a 
single  peak  compressive  load  followed  by 
several  smaller  peaks.  For  theses  cases,  any 
hersteresis  in  cyclic  loading  after  the  initial 
pulse  is  generally  viewed  as  unimportant. 

In  vibration  induced  sinusoidal  motion  where 
cyclic  shear  is  a  predominant  effect,  hysteresis 
behavior  caused  by  cyclic  loading  must  be 
considered.  In  order  to  model  this  energy 
dissipation  process,  a  viscoelastic  stress-strain 
relation  is  incorporated  into  the  CAP  model  to 
simulate  the  effects  of  low  amplitude  cyclic 
hysteresis  within  the  yield  surface.  This  allows 
for  a  controlled  amount  of  energy  dissipation  in 
cyclic  loading  and  can  be  fine  tuned  to  provide 
the  correct  damping  at  the  frequencies  asso¬ 
ciated  with  the  forcing  functions. 

Experimental  results  indicate  that  the  amount 
of  hysteresis  for  soils  is  relatively  rate  indepen¬ 
dent.  Because  of  this  experimental  observation, 
the  following  approach,  as  discussed  in  detail  in 
Reference  (9)  is  taken; 


Linear  viscous  damping  is  included  to  model 
stress  paths  within  the  yield  surface  of  the  CAP 
model.  A  standard  solid  is  chosen  to  represent 
the  shear  behavior  for  stress  paths  within  the 
current  yield  surface.  This  standard  solid  model 
requires  three  material  parameters  to  describe 
the  deviatoric  viscoelastic  response.  These 
parameters  are  an  instantaneous  modulus  Gf,  a 
long  term  modulus  Gs  and  a  relaxation  frequen¬ 
cy  to. 

The  viscoelastic  portion  of  the  model  is  de¬ 
scribed  by 

S  +  cos  =  IGf^  +  'bDGsf  (I) 


where  s  and  eE  are  the  stress  and  elastic  strain 
deviators,  is  a  relaxation  frequency  and  Gf  and 
Gs  represent  the  shear  moduli  under  fast  and 
slow  loading  respectively. 

The  volumetric  behavior  is  elastic 

Jl  =  3  (2) 


£• 

where  J  i  and  e  v  are  the  traces  of  the  stress  and 
strain  tensors  and  K  is  the  bulk  modulus. 

To  obtain  numerical  values  for  these  parame¬ 
ters  an  experimentally  justified  approximation 
is  made  that  both  the  hysteresis  and  the  average 
modului  are  relatively  frequency  independent. 
The  two  parameters  G,  and  Gf  are  related  to  an 
average  modulus  over  a  loading  cycle  and  to  the 
energy  loss  per  cycle  at  a  central  frequency  of 
interest.  The  constants  Gs  and  Gf  can  be 
obtained  from  a  single  load  -  unload  curve  of 
the  compacted  native  backfill.  The  parameters 
Gs  and  Gf  are  related  to  the  average  modulus 
over  the  cycle  and  to  the  energy  loss  per  cycle. 
The  parameter  to  is  introduced,  where  is  the 
circular  frequency  at  which  the  energy  loss  per 
cycle  of  a  standard  solid  reaches  a  maximum. 
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The  parameter  o>  is  chosen  to  correspond  to  the 
central  frequency  of  interest,  that  is  the  forcing 
function  frequency.  Thus  for  the  main  frequen¬ 
cy  of  interest,  both  the  hysteresis  and  the 
average  modului  are  relatively  frequency  inde¬ 
pendent. 

The  viscoelastic  algorithm  begins  by  comput¬ 
ing  a  viscoelastic  “trial  value”  for  the  stresses, 
i.e.,  a  value  based  on  the  assumption  that  no 
failure  or  cap  plasticity  occurs.  If  the  assump¬ 
tion  is  found  to  be  false,  the  stresses  are  revised 
by  means  of  a  “plastic  correction”  (Ref.  10). 

For  reasons  of  accuracy  under  arbitrary  loading 
paths,  the  viscoelastic  trial  values  are  obtained 
for  a  strain  increment  Ae  during  a  time  step  At 
as  the  exact  solution  for  linear  straining  during 
the  time  step  t„  <  t  ^  t„  +  At, 


at  t  =  t„  +  At: 


5"+1  =  Sne~a>Lt  +  2Gs*£"(1  -  e-"*) 

+  2[GS  +  (i  _  e-wi“)WE 
w  At 

(6) 

which  is  the  viscoelastic  trial  deviatoric 
stresses. 

The  viscoelastic  trial  stresses  are  then  checked 
against  the  failure  envelope  and  cap.  If  either  or 
both  of  these  surfaces  is  exceeded,  the  stress 
point  is  corrected  so  as  to  satisfy  the  appropriate 
yield  condition  and  normal  flow  rule. 


<*«>-«*■  +  *. <3> 

where  eE(t)  is  the  viscoelastic  deviatoric  strain 
at  tN-  Substituting  equation  3  into  equation  1 
produces 

A<?  „  „  -  A(?  / 

s  +  tos  =  2 Gf  — — +  2mGs  [eE’'  +  — — {t- ^)j 
At  At 

(4) 

which  may  be  solved  using  s(tn)  =  sn,  to  give 


The  viscoelastic  model  was  fit  to  the  available 
data  for  the  backfill  soil  (Figure  (3). 

The  bulk  modulus  K  was  estimated  from  the 
uniaxial  stress-strain  properties  of  the  com¬ 
pacted  native  backfill. 

The  “fast”  and  “slow”  shear  moduli,  Gf  and  Gs 
respectively  were  related  to  the  size  of  the 
elliptic  hysteresis  loop  and  are  dependent  upon 
the  amount  of  energy  dissipated  through  vis¬ 
cous  damping  in  the  soil  material. 


The  viscoelastic  relaxation  parameter  oo  is 
taken  to  be  equal  to  the  forced  frequency  of 
excitation  of  the  structure. 

For  the  backfill  soil  material,  the  properties 
used  to  fit  the  viscoelastic  CAP  model  are 
estimated  from  the  uniaxial  stress-strain  prop¬ 
erties  of  the  compacted  native  backfill  (Fig.  3) 
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exhibiting  a  confined  modulus  (K  +  4/3  G)  of 
26000  psi  and  are  as  follows: 

K  =  18000.  psi 

Gs  =  6000.  psi 

Gf  =  9000.  psi 

CW  =  0.05 

CR  =  3.0 

CA=  1.02  ksi 


CB  =  0.2  (1/ksi) 
CC  =  1 .0  ksi 
CD  =  1.5  (1/ksi) 
EL  =  0.0  ksi 
LTYPE  =  1 
TCUT  =  .02  ksi 
FCUT  =  .02  ksi 
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SECTION  3 


RESULTS  OF  BURIED  ARCH  SIMULATION  STUDIES 


3.1  DISCUSSION  OF  RESULTS. 

The  physical  results  obtained  from  the  ESSEX 
V  buried  structure  tests  differed  appreciably  in 
form  from  experimental  predictions.  In  an 
attempt  to  explain  the  reasons  for  this  differ¬ 
ence,  various  damping  mechanisms  have  been 
identified  as  possible  causes  of  discrepancy 
between  experimental  and  analytical  results 
(Ref.  2).  These  damping  mechanisms  include: 

1.  Radiation  damping 

2.  Interface  friction  between  the  soil  and  the 
structure 

3.  Hysteresis  in  the  soil  backfill  adjacent  to 
the  structure 

4.  Structural  damping 

Previous  work  (Ref.  2)  focused  on  the  mecha¬ 
nisms  of  radiation  damping  and  friction,  as  they 
were  previously  believed  to  be  capable  of  being 
the  major  contributors  to  the  total  effective 
damping. 

The  prese  it  study  focuses  on  the  role  of  soil 
hysteresis  in  contributing  to  the  total  effective 
damping  needed  to  supress  resonance.  Soil 
hysteresis  causes  energy  to  be  absorbed  during 
the  cycles  of  load-unload-reload  produced  by 
the  driving  force  oscillations.  The  actual 
amount  of  effective  hysteresis  damping  is  a 
function  of  the  amount  of  energy  absorbed  in 
each  cycle. 

Mechanical  impedance  data  obtained  from  the 
ESSEX  V  test  series  showed  the  variation  in 
impedance  levels  (F/V)  due  to  changes  in  the 
forcing  function  frequency.  Each  trough  in  the 
impedance  curve  represents  a  frequency  at 


which  resonance  occurs  and  represents  a  natu¬ 
ral  frequency  of  the  buried  structure. 

Analytical  impedance  results  were  obtained  for 
the  three  test  series  which  used  a  500  pound 
(2224  N)  forcing  function.  These  series  are; 

1.  A  single  vibrator  at  the  arch’s  crown 

2.  Two  vibrators  in  phase 

3.  Two  vibrators  out  of  phase 

For  each  test  series,  a  frequency  sweep  of  50  to 
1000  hertz  was  experimentally  performed. 

The  analytical  simulation  studies  were  per¬ 
formed  for  a  frequency  sweep  from  50  to  250 
hertz  in  intervals  of  10  hertz.  For  each  frequen¬ 
cy  in  the  interval,  a  direct  integration  in  time 
was  performed.  That  is,  21  separate  time 
integration  analyses  were  performed  to  obtain 
the  frequency  sweep.  A  more  limited  frequency 
range  was  used  in  the  analytical  studies  versus 
the  experimental  studies  because  of  limitations 
in  the  analytical  procedure  at  both  low  (less 
than  50  hertz)  and  high  (greater  than  250  hertz) 
frequencies.  The  reasons  for  this  are  that  in  the 
low  frequency  range,  the  use  of  absorbing 
boundaries  along  the  finite  element  mesh 
boundaries  restricts  the  rigid  body  motion 
prevalent  for  these  low  frequencies,  thereby 
invalidating  the  analytical  procedure  and  for 
the  high  frequency  range,  the  size  of  the  finite 
element  mesh  restricts  the  frequency  content  of 
the  solution,  thereby  invalidating  the  analytical 
procedure. 

The  simulation  study  considers  a  cross-section 
of  the  arch  with  length  effects  assumed  to 
satisfy  plane  strain  kinematics.  The  forcing 
function  for  the  plane  strain  analyses  is  force 
per  unit  length  of  the  arch  (F/L).  Therefore,  the 
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analytically  obtained  impedance  values  differ 
from  those  measured.  The  analytical  imped¬ 
ance  values  are  obtained  from  (F/L)/V.  The 
analytical  impedance  values  differ  in  units 
from  the  experimental  values  by  a  inverse 
length  unit.  This  discrepancy  can  be  eliminated 
by  assuming  a  characteristic  length  L*  for 
which  the  plane  strain  analyses  impedance  can 
compare  with  the  full  three  dimensional  re¬ 
sponse.  A  possible  value  for  L*  could  be  the 
arch’s  radius.  Then  the  analytical  impedance 
values  would  be  obtained  from  (F/V)*(L*/L). 

To  isolate  the  effects  of  soil  hysteretic  nonlin¬ 
earity  upon  impedance  versus  forcing  function 
frequency  relations,  results  are  presented  for 
two  assumed  backfill  soil  models.  The  first  soil 
model  is  an  elastic -plastic  CAP  model  and  the 
second  is  a  viscoelastic  CAP  model. 

For  the  single  vibration  test,  using  a  500  lb 
(2224  N)  force  at  drive  point  FI  (crown  of 
arch),  the  effect  of  backfill  soil  hysteresis  is 
illustrated  in  figure  (4).  For  the  case  of  elas¬ 
tic-plastic  backfill  soil  model,  two  distinct 
troughs  in  the  impedance  curve  are  obtained. 
These  occur  at  100  and  180  hertz  and  represent 
frequencies  where  structural  resonance  occurs. 
However,  the  use  of  a  viscoelastic-plastic 
backfill  model  supresses  the  resonance  occur- 
ing  at  1 80  hertz  and  smooths  out  the  impedance 
trough  which  occurs  at  100  hertz.  Thus  the  use 
of  a  viscoelasticity  backfill  CAP  model  reduces 
the  tendency  of  the  structural  arch  to  resonate 
within  the  frequency  band  considered. 

The  comparisons  of  both  analytical  results  with 
the  experimental  data  for  the  single  vibration 
test  are  shown  in  figures  5  and  6.  The  exper¬ 
imental  response  of  the  two  gages,  CL6  and 
CL8  are  shown.  The  response  of  the  gages  CL6 
and  CL8  illustrate  the  change  in  impedance 
along  the  length  of  the  arch’s  crown.  If  the  arch 
responded  in  plane  strain  modes,  these  two 
gages  should  have  identical  responses.  This  is 
not  observed  as  gage  CL8  has  a  greater 


impedance  than  gage  CL6  because  it  is  further 
removed  from  the  force  drive  point.  For  the 
elastic  backfill  soil  model,  the  computed  reso¬ 
nance  at  180  hertz  is  not  observed  in  the 
experiment  at  either  of  the  measured  gages.  For 
the  viscoelastic  (hysteretic)  backfill  soil  model, 
the  computed  response  shows  a  resonance  at  90 
hertz.  A  resonance  of  80  hertz  is  exhibited  by 
gage  CL8.  Note  that  the  magnitude  of  the 
calculated  response  differs  from  the  magnitude 
of  the  measured  response  because  of  the 
analytical  assumption  of  plane  strain  response. 
What  is  critical,  however,  is  that  the  inclusion 
of  a  hysteretic  backfill  soil  model  alters  the 
shape  of  the  response  impedance  diagram  by 
supressing  the  resonance  mode  at  180  hertz. 

For  the  dual  in-phase  vibration  test,  using  a  500 
lb  (2224  N)  force  at  drive  points  F2  and  F3,  the 
effect  of  backfill  soil  hysteresis  is  illustrated  in 
Figure  (7).  For  the  case  of  elastic-plastic 
backfill  soil,  a  distinct  trough  in  the  impedance 
curve  is  obtained.  This  occurs  at  140  hertz  and 
represents  structural  resonance.  The  use  of  a 
viscoelastic-plastic  backfill  soil  model  does  not 
supress  the  resonance  occuring  at  140  hertz. 

The  comparisons  of  both  analytical  results  with 
the  experimental  data  for  the  dual  in-phase 
vibrator  test  are  shown  in  Figures  (8)  and  (9). 
The  experimental  response  of  the  gage  A3  is 
shown.  Note  that  the  magnitude  of  the  com¬ 
puted  impedance  curve  differs  from  the  magni¬ 
tude  of  the  measured  response  because  of  the 
plane  strain  limitation.  As  seen  in  Figure  (7), 
the  inclusion  of  the  hysteretic  soil  model  did  not 
supress  the  calculated  resonance  mode  at  140 
hertz.  This  resonance  mode  was  not  observed  in 
the  experimental  records. 

For  the  dual  out-of-phase  vibrator  test,  using  a 
500  lb  (2224  N)  force  at  drive  points  F2  and  F3, 
the  effect  of  backfill  soil  hysteresis  is  illustrated 
in  Figure  (10).  For  the  case  of  elastic-plastic 
backfill  soil,  a  distinct  trough  in  the  impedance 
curve  is  obtained.  This  occurs  at  1 10  hertz  and 
represents  a  frequency  where  structural  reso¬ 
nance  occurs.  However,  the  use  of  a  viscoelas- 
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tic-plastic  backfill  soil  model  suppresses  the 
resonance  occuring  at  110  hertz.  Thus  the  use 
of  a  viscoelastic  backfill  CAP  model  reduces 
the  tendency  of  the  structural  arch  to  resonate 
within  the  frequency  band  considered. 

The  comparisons  of  both  analytical  results  with 
the  experimental  data  for  the  dual  out-of-phase 
vibrator  test  are  shown  in  Figures  (11)  and  (12). 
The  experimental  response  of  gage  A3  is 
shown.  For  the  elastic  backfill  soil  model,  the 
computed  resonance  at  110  hertz  is  not  ob¬ 
served  in  the  experiment  gage  record.  For  the 
viscoelastic  (hysteretic)  backfill  soil  model,  the 
computed  response  supresses  the  resonance  at 
110  hertz.  The  inclusion  of  a  hysteretic  backfill 
soil  model  alters  the  shape  of  the  response 
impedance  diagram  by  supressing  the  reso¬ 
nance  mode  at  110  hertz. 


For  the  single  vibrator  test  and  the  dual 
out-of-phase  vibrator  test,  the  inclusion  of 
viscoelasticity  into  the  backfill  soil  model 
served  to  supress  resonant  modes.  The  supres- 
sion  of  these  resonance  modes  agreed  with  the 
experimental  data  which  showed  that  for  the 
arch  structure,  resonance  was  supressed  in  the 
covered  tests  as  compared  to  the  uncovered 
tests. 

For  the  dual  in-phase  vibratortest,  the  inclusion 
of  viscoelasticity  into  the  backfill  soil  model 
did  not  serve  to  supress  the  resonant  modes. 
This  supression  was  observed  however  in  the 
experimental  gage  records.  Possibly,  the  reason 
for  this  difference  is  that  the  plane  strain 
assumptions  were  not  at  all  valid  for  this 
in-phase  test  series. 
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SECTION  4 


CONCLUSIONS  AND  RECOMMENDATIONS 


4.1  CONCLUSIONS. 

1.  Resonant  frequencies  for  the  covered 
arch  with  a  single  vibrator  force  at  drive 
point  F 1  are  analytically  supressed  when 
a  viscoelastic-plastic  soil  backfill  CAP 
model  is  used. 

2.  Resonant  frequencies  for  the  covered 
arch  with  a  dual  out-of-phase  vibrator 
force  at  drive  points  F2  and  F3  are 
analytically  supressed  when  a  viscoelas¬ 
tic-plastic  soil  backfill  CAP  model  is 
used. 

3.  Resonant  frequencies  for  the  covered 
arch  with  a  dual  in-phase  vibrator  test  at 
drive  points  F2  and  F3  are  not  analytical¬ 
ly  supressed  when  a  viscoelastic-plastic 
soil  backfill  CAP  model  is  used. 

4.  The  inclusion  of  a  viscoelastic-plastic 
soil  model  is  critical  in  helping  to 
explain  the  resonant  mode  supression 


observed  in  the  covered  arch  experimen¬ 
tal  results. 

5.  Impedance  testing  excites  nonlinear  soil 
response  typified  by  hysteretic  energy 
losses  and  radiation  and  friction  energy 
dissipation. 

6.  The  applicability  of  impedance  testing  to 
define  buried  arch  response  is  question¬ 
able. 

4.2  RECOMMENDATIONS. 

Additional  experimental  and  analytical  studies 
are  needed  to  expand  the  data  base  of  buried 
structural  response.  Soil  parameters  defining 
hysteretic  energy  losses  need  to  be  more 
accurately  defined  to  assist  in  establishing 
viscoelastc-plastic  soil  constitutive  models  to 
simulate  the  hysteretic  energy  losses  which 
occur  during  the  forced  vibration  tests  of  buried 
structures. 
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Figure  2.  Finite  element  model  of  buried  arch. 
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Figure  3.  Uniaxial  stress-strain  properties  of  compacted 
native  backfill  ESSEX  V  arch. 
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Figure  4.  ARCH  structure  3A,  analytical  response  at  crown  -  single  vibration  test 
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Figure  5.  Elastic  soil  cap  model  -  ARCH  structure  3At 

analytical  response  at  crown  -  single  vibration  test, 
FI  =  500  lbs  (2224  N). 
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Hysteretic  soil  cap  model  -  ARCH  structure  3A, 
analytical  response  at  crown  -  single  vibration  test, 
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Figure  7.  ARCH  structure  3A,  analytical  response  at  gage  location  A3, 
single  vibration  test,  FI  =  500  lbs  (2224  N). 
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Figure  9.  Hysteretic  soil  cap  model  -  ARCH  structure  3A, 

response  data  at  gage  A3  -  two  vibrators  in  phase. 
F2  =  F3  =  500  lbs  (2224  N). 
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NOTATION 


CA,  CB  CC,  CD,  CW,  FCUT,  TCUT  CAP  fit  Parameters 
e®  elastic  strain  devisor 

trace  of  strain  tensor 

fMAX  maximum  harmonically  applied  force 
Gf  shear  modulus  under  fast  loading 

G,  shear  modulus  under  slow  loading 

Ji  trace  of  stress  tensor 

K  bulk  modulus 

S  stress  deviator 

VMAX  maximum  velocity 

(0  relaxation  frequency 
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